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Abstract 

The quasi-parallel photon-photon scattering by combining two-color laser fields is an approach to 
produce resonant states of low-mass fields in laboratory. In this system resonances can be probed 
via the four-wave mixing process in the vacuum. A search for scalar and pseudoscalar fields was 
performed by combining 9.3 /uJ/0.9 ps Ti-Sapphire laser and 100 /rJ/9 ns Nd:YAG laser. No 
significant signal of four-wave mixing was observed. We provide the upper limits on the coupling- 
mass relation for scalar and pseudoscalar fields, respectively, at a 95% confidence level in the mass 
region below 0.15 eV. 
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FIG. 1: Quasi-parallel colliding system by combining two-color laser fields where beam diameter 
d, focal length /, and the incident angle takes 0 < ?? < A9 which is unavoidable due to ambiguity 
of the wave vectors of incident photons by the nature of focused lasers. 


INTRODUCTION 


Uncovering the nature of dark energy and dark matter is one of the most crucial problems 
in modern physics. Low-mass and weakly coupling helds predicted by theoretical models in 
cosmology and particle physics can be candidates for such dark components. For instance, 
based on the scalar-tensor theory with the cosmological constant A (STTA) [l|, dark en¬ 
ergy is interpreted as decaying A while the universe becomes older due to the gravitational 
coupling between extremely light dilatons, a kind of scalar helds (</>), and matter helds. Ob¬ 
serving the 77 ^ 0 —>■ 77 process with extremely high intensity laser helds can be a method 
of searching for 0 in laboratory j^. The same approach can also be applied to searches for 
ow-mass pseudoscalar helds (a), if the photon spin states are properly chosen [^. Axion 
4, 5], a pseudoscalar held associated with breaking of Peccei-Quinn symmetry j^, is a suit- 


Dle. Axion is supposed to be one 
7|, l8(]. Therefore, these theoretical 


able candidate to which this method is directly applica 
of the most reasonable candidates for cold dark matter 
models strongly motivate us to search for such helds in laboratory in general. 

Axion searches via the two photon coupling processes have been performed by a number 


of experiments, for example, solar axion searches 9l-ll5l|. light shining through a wall 16l-ll9|. 


and the axion dark matter experiment 291. l2ll| . Following the hrst search for scalar helds 
at quasi-parallel colliding system (QPS) j^, the upgraded search for sub-eV scalar and 
pseudoscalar helds is presented in this paper. 
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With the schematic view of QPS in Fig. [H we briefly explain the essence of onr method 
as follows. By using variables dehned at QPS, the center of mass system (CMS) energy 
between a randomly selected photon pair is expressed as 

^CM 5 = 2a; sin 7?, (1) 


where a; is the energy of incident photons and'd is half of the incident angle of the photon 
pair. Extremely low collision energies are realizable at QPS by focusing a laser held because 
small values of '& can be automatically introduced. 

In order to overcome low scattering amplitudes of 77 ^ 0 /ct —>■ 77 processes due to weak 
coupling, we hrst utilize the character of the integrated resonance effect by capturing Ecms 
within AEcms via A9 prepared by a creation laser held. Secondly we let another laser 
held propagate into the optical axis common to the creation laser. This laser induces decay 
of resonance states into a specihc energy-momentum space by the coherent nature of the 
inducing held. The scattering probabil^ is thus proportionally increased by the number of 


photons in the inducing laser held j2 , l3|, l23|, l24 1 


Energies of decayed photons are dehned by the following energy conservation 


a; -1- a; = (2 — u)oj + uoo, (2) 

where u is an arbitrary number which satishes 0 < m < 1. We re-dehne the energies of hnal 
state photons as following 


a;3 = (2 — u)oj, 

Ui = no;, (3) 

where a ;3 and a ;4 are energies of signal photon and inducing photons, respectively. 

In the case of the scalar held exchange, the relation of linear polarization states between 
initial and hnal state photons when the wave vectors are on the same reaction plane are 
expressed as follows: 

a;{l}-|-a;{l} = a; 3 {l}-|-a; 4 {l}, 

a;{l}-|-a;{l} = a;3{2} -|- a;4{2}, (4) 

where {1} and {2} are linear polarization states orthogonal to each other. In the pseu¬ 
doscalar hied exchange, the polarization relation are expressed as 
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a;{l} + cu{2} — cjajl} + U4{2}, 
Cj{l} + Cj{2} = CJ3{2} + CJ4{1}. 


( 5 ) 


We emphasize that above relations are limited only to the theoretically ideal case where all 
fonr photons are on the same reaction plane within the treatment based on plane waves. In 
the focnsed QPS, however, we must accept independent rotations of the incident pi—p 2 plane 
and the outgoing p^, — p^ plane as illustrated in Fig. [13] with respect to an experimentally 
given linear polarization plane. This implies that even if we supply u as the pure {l}-state 
by a polarizer at the moment of plane wave propagation in advance of focusing, mixing of 
{1} and {2} states for randomly selected incident photon pairs is unavoidable while lasers 
are focused. Therefore, the focused QPS with a hxed initial linear polarization plane has 
sensitivity to both scalar and pseudoscalar helds simultaneously. We discuss about this 
nature in detail in Appendix A. 

The relation in Eq.(l2]) is similar to ’’four-wave mixing” in matter corresponding to the third 
order non-linear quantum optical process in atoms 25|,l26|]. Therefore, the observation of the 
four-wave mixing process in the vacuum may be interpreted as a replacement of the atomic 
nonlinear process by the exchange of unknown scalar or pseudoscalar helds. The observation 
of four-wave mixing in the vacuum is also used as a method for testing higher-order QED 


effect 


27 


3Q|. 


Photons produced via the atomic four-wave mixing process can be the main background 
source for this search. The hrst search for scalar helds at QPS 22| was performed with weak 
intensity lasers, thus, the ehect of the four-wave mixing process in atoms was negligible. In 
this experiment, however, the four-wave mixing photons originating from the residual gas 
are anticipated due to much higher beam intensities. In this paper the method to obtain 
the exclusion limits in the search at QPS sensitive to both scalar and pseudoscalar helds 
is provided under the circumstance where a hnite amount of background photons must be 
evaluated. 
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II. THE COUPLING-MASS RELATION 


The effective interaction Lagrangians coupling between two photons and 0 / cr are expressed 
as 

-L^ = -L, = gM-^^F,,F^’^a, ( 6 ) 

where M has the dimension of energy and g is the dimensionless constant. The yield of 
signal photons, y, is expressed with experimental parameters relevant to lasers and optical 
elements as follows: 


3^ = 


A. 


64\/27r'^ \ CTc 


5 14A 


Tid? \ (u — uy ( gm[eV] A ^ /m[eV] \ ^ 


uu \M[eV]J V^[eV] 


WOFsC^bN^N,, 

( 7 ) 


where the subscripts c and i indicate creation and inducing laser, respectively, A is wave¬ 
length, r is pulse duration, / is focal length, d is beam diameter, u and u are upper and 
lower values on u determined by the spectrum width of respectively, m is mass of the 
exchanging held, W is the numerical factor relevant to the integral of the weighted resonance 
function which is rehned in Eq. flSOp in Appendix B compared to W ~ 7r/2 in Ref. 2^, Q 
is the incident-plane-rotation factor described in Appendix A, Fs is the polarization de¬ 
pendent axially asymmetric factor for outgoing photons [s], Cmh is the combinatorial factor 
originating from selecting a pair of photons among multimode frequency states and N is 
the average numbers of photons in the coherent state. The detail of the formulation of the 
signal yield is summarized in Appendix of Ref. 2^. The coupling constant g/M is expressed 


as 


M[eV] 


= 2^/^87r“ 




>>w3[eV] 


(i) te) (h (i7&) 


(r— 


WGFsCmbNyN, 


-m ^/^[eV]. (8) 


III. EXPERIMENTAL SETUP 


We explain the experimental setup to detect signals of four-wave mixing in the vacuum. 
The schematic view of the setup is shown in Fig. O 

A Ti-Sapphire laser (wavelength 800 nm) and a Nd:YAG laser (wavelength 1064 nm) are 
used as the creation and the inducing lasers, respectively. To reduce the number of back- 
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ground photons emitted from the residual gas via four-wave mixing, the linear polarization 
states of the creation and inducing lasers are conhgured to linear polarization states {1} 
and {2}, respectively. The beam alignments of the lasers are monitored by CCD cameras 
(CCD) and the pulse energies of the creation and inducing lasers are measured by photo¬ 
diodes (PD). These beams are combined by a dichroic mirror (DM). The combined beams 
are guided into the vacuum chamber at the 20 mm beam diameter and focused with the 
convex lens at the focal length 200 mm. 



FIG. 2: The schematic view of the experimental setup 

The expected wavelength of the corresponding signal photon is evaluated from the following 
equation 

A, - XJ2 ^ ^ 

A light source with the central wavelength of 640 nm is combined with the creation and 

inducing lasers by DM to evaluate the detection efficiency and to trace the trajectory of 
signal photons for the detector alignment. 

The agreement of the optical axes between the two lasers are adjusted at a precision of 2-3 
fim by monitoring individual beam prohles at the near side and the far side of the focal spot 
with the CCD camera. The beam prohles at the focal spot are shown in Fig. |3l The spot 
sizes of the creation and inducing lasers which are dehned as 2 a of the 2D Gauss functions 
htting the beam prohles, are 21 pm and 23 pm, respectively. The creation laser overlaps 
with 87 % of the beam energy of the inducing laser at the focal spot. Thus, the ehective 
beam energy of the inducing laser is evaluated by correcting the measured beam energy with 
this overlapping factor. 
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FIG. 3: The beam profiles of the creation laser (left) and the inducing laser (right) at the common 
focal point captured by a common CCD camera. 


Signal photons generated within the focal volnme travel along the common optical axis of 
the combined lasers. Signal photons are separated from the creation and indncing lasers by 
the prism and signal wave hlters are placed to farther eliminate the residnal photons from 
the combined lasers. The polarization beam splitter (PBS) transmits {l}-polarized photons 
and reflects {2}-polarized photons. Incident photons are split between the shorter optical 
hber Path{l} and the longer Path{2}. The incident photons to PBS are eventually observed 
by the common photo-device having relative time delay of 23 ns. We use a single-photon- 
countable photomultiplier tube (PMT) R7400-01 manufactured by HAMAMATSU as the 
photo-device. 

The repetition rate of the creation laser is 1 kHz and that of the inducing laser is 10 Hz by 
synchronizing the trigger with the 1 kHz creation pulsing. The data acquisition trigger of 
20 Hz is synchronized with the 1 kHz creation laser pulsing which includes pedestal triggers 
in order to provide four patterns of triggers. The time coincidence between creation and 
inducing pulses are performed by adjusting the relative injection timing between the two 
lasers so that the relative time maximizes the four-wave mixing yield in the air. The shutter 
is placed on the creation laser beam line and it repeats open and close every 5 sec. We 
acquire data with the four patterns of triggers, which are ’’both of lasers are incident (S)”, 
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’’only creation laser is incident (C)”, ’’only inducing laser is incident (I)”, and ’’neither of 
lasers are incident (P)”. The digital oscilloscope recorded waveform data from the PMT 
and two photo-diodes synchronized with the 20 Hz data acquisition trigger. The recorded 
waveform data from the PMT are sorted into four types of trigger patterns S, C, I and 
P. The four trigger patterns are classihed by checking the charge correlations between the 
waveform data from the two photo-diodes for intensity monitoring. 


IV. METHOD OF THE WAVEFORM ANALYSIS 

The observed photon counts are estimated by analyzing the waveform data from the PMT. 
The individual waveform consists of 500 sampling data points within a 200 ns time window. 
We search for negative peaks of which amplitude exceed a given threshold. We then calculate 
charge sums of the peak structures. Figure 0] shows a sample of waveform data where peak 
structures are identihed. Charge sums of peak structures are evaluated in units of the 
single-photon equivalent charge, —4.21 x 10“^^C. 



FIG. 4: The waveform data sample which has two peak structures. Black shaded areas show the 
integral ranges to evaluate charge sums of individual peak structures. 

There are some accidental noisy events among recorded waveform data. In our analysis 
method, these noise structures could be misidentihed as large photon-like peak structures. 
Therefore, it is necessary to remove such noisy events from analyzed waveform data before 
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counting photon-like peaks. We can identify noisy events by analyzing the frequencies of 
waveforms. Noisy waveforms tend to have lower frequencies than those of normal waveforms. 
The frequencies are estimated by counting the number of nodes which is defined as the 
intersections between a waveform and the average line of amplitudes within the 200 ns time 
window. The distributions of the number of nodes for each trigger pattern are shown in 
Fig. |5l We regard a waveform of which the number of nodes is lower than 150 as a noisy 
event in all trigger patterns by conhrming that the differences of the distributions among 
four trigger patterns are not prominent. The typical waveforms of noisy events and normal 
events identihed by this method are shown in Fig. [HI 

s c 






FIG. 5: Distributions of the number of nodes for trigger patterns S,C,I and P. The events with the 
fewer number of nodes below the red vertical line are identified as noisy events. 


V. MEASUREMENT OF THE FOUR-WAVE MIXING PROCESS IN THE 
RESIDUAL GAS 

The background photons can be produced via the four-wave mixing process occurred in 
residual atoms in the vacuum chamber. To estimate the expected number of background 
photons, we measured the pressure dependence of the number of four-wave mixing photons 
in gas. Figure [7] shows arrival time distributions of observed photons in the air at 5.0 x 
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FIG. 6: Examples of waveforms of noisy events and normal events. The two panels in the left 
and right sides show waveforms of normal events and noisy events, respectively. The red horizontal 
lines indicate the averages of the amplitudes of sampling points for each waveform. 


Pa among four trigger patterns. Specific two peak structures appear only at S-pattern. 
These peak structures have approximately 23ns time interval, which agrees with the optical 
path-length difference between Path{l} and Path{2}. We count the number of photons 
within a time domain T{1} (71-75 ns) with {l}-polarized state and T{2} (94-98 ns) with 
{2}-polarized state. 

The number of four-wave mixing signals Ns are evaluated from the following equation. ( 
see Eqs.(18) and (19) in Ref. 22|) 


.. hP5 ^ VPs 

Ns = ns- r—nc - —ni + —rip, 

Wc Wi Wp 


( 10 ) 


where n* and Wi denote the number of photon-like peaks in the signal domains and the 
number of events in trigger pattern i, respectively. 

The pressure dependence of the number of four-wave mixing photons per S-trigger event are 
shown in Fig. [HI Data points are fit by the quadratic function of pressure. We extrapolate 
the number of four-wave mixing photons in the residual gas at 2.3 x 10“^ Pa (an equivalent 
condition to the vacuum data we discuss later) from the fitting function. The efficiency- 
corrected number of {l}-polarized and {2}-polarized photons in residual gas Afgasi and Mgas 2 
with the same shot statistics as the vacuum data are evaluated as follows: 


10 














s 




c 




FIG. 7: The arrival time distributions of observed photons per triggered event (efficiency- 
uncorrected) at 5.0 x Pa. The left and right bands bounded by neighboring two red lines 
in each panel indicate the time domains T{1} and T{2} where {1} and {2}-polarized photons are 
expected to be observed, respectively. In this figure, the threshold value for peak identification are 
set lower than that of the actual data analysis on purpose to show typical pedestal structures in 
each trigger pattern. 


Mgasl = 1.7 ±1.1 X 10-5, 

= 1.7 ±1.1 X 10-5. (11) 

We confirmed that the expected value of four-wave mixing photons from the residual gas 
are negligibly small in the vacuum data for a given total statistics. 

VI. SEARCH FOR FOUR-WAVE MIXING SIGNALS IN THE VACUUM 

We acquired data at 2.3 x 10-^ Pa for the search for the resonant states of 0 and a helds. 
Figure |9] shows the arrival time distributions of observed photon counts . Table [I] summarizes 
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FIG. 8: The pressure dependence of the number of four-wave mixing photons in the residual gas 
inside the interaction chamber per S-trigger event. The red and black lines represent the htting 
functions for {1} and {2}-polarized states, respectively. 

the numbers of observed photon-like signals evaluated in units of the single-photon equivalent 
charge with {1} and {2}-polarized states for each trigger pattern, respectively. 

TABLE I: The numbers of observed photons in Tj!} and r{2} for each trigger pattern, nn and nj 2 
are the number of photons evaluated in units of single-photon equivalent charge in trigger pattern 
i with {1} and {2}-polarized states, respectively. Wi is the number of events in trigger pattern i. 


Trigger i nn 

ni2 Wi 

S 

0 

0 46120 

C 

0 

0 46203 

I 

0 

0.07 46044 

P 

0 

1.53 46169 


After performing subtractions between four patterns of histograms in Fig. |9] based on the 
relation in Eq. ffTOl) . we obtained the time distribution of Ns as shown in Fig. [101 The 
number of signals with {1} and {2}-polarized states are, respectively, given as follows: 
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FIG. 9: Arrival time distributions of observed photons at 2.3 x 10 ^ Pa. The data points in each 
trigger pattern are normalized to the number of triggered events of S-trigger pattern. 


Nsi = 0 ± 0(stat.) ± 2.16(syst.I) ± 0.30(syst.n) ± O(syst.in), 

Ns 2 = 1-46 ± 1.27(stat.) ± 2.16(syst.I) ± 0.04(syst.n) ± 3.59(syst.in). (12) 

The systematic error I originates from the number of the photons out side of the two arrival 
time windows for {1} and {2}-polarized states. This was evaluated by calculating the root 
mean square of Ns except in the T{1} and T{2} windows. The systematic error If originates 
from the dependence on the threshold values for the peak finding —1.3 ± 0.1 mV. The 
systematic error hi is relevant to the ambiguities of the rejection of noisy events 150±5 nodes. 


VII. THE EXCLUDED COUPLING-MASS LIMITS FOR SCALAR AND PSEU¬ 
DOSCALAR FIELDS 

There is no signihcant four-wave mixing signal in this search from the result in ffl^ . We 
thus evaluate the exclusion regions on the coupling-mass relation as follows. 
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FIG. 10: The arrival time distribution of Ns defined in Ea. (|10p . 


We estimate the upper limit on the sensitive mass range as 


m < 2ci;sinA6* 



0.15 eV 


(13) 


based on values summarized in Table UTl where'd in Fig. [T] varies from 0 to A6 dehned by a 
focal length / and a beam diameter d. 

The number of efficiency-corrected {l}-polarized signal photons Afsi and that of {2}- 
polarized signal photons A /52 are evaluated from the following relations with the experi¬ 
mental parameters 


Ms, = Ms, = (14) 

^optl^D ^opt2^D 

where eopti and eopt 2 are the attenuation ratios of the signal photons propagating from the 
interaction point through Path{l} and Path{2}, respectively. 

These attenuation factors are composed of the transmittance of optical devices and the ac¬ 
ceptance of signal paths with respect to the actual location of the PMT. They are inclusively 
evaluated by sampling the beam energies of the 640 nm calibration light at the focal point 
and the detection point, respectively, and taking the ratio between them. The matching 
of beam paths between the calibration light and four-wave mixing signals are ensured by 
adjusting the beam center of calibration light with respect to those of creation and induc¬ 
ing lasers at the near side and the far side of the focal spot, respectively, e^i is the signal 
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TABLE II: Data table of experimental parameters. Qff and Q 12 represent the incident-plane- 
rotation factor for scalar and psendoscalar field exchanges, respectively. The evaluation of Q is 
discussed in Appendix A. ^^1122 ^-nd ^^1212 denote the axially asymmetric fractor for scalar and 
psendoscalar field exchanges, respectively. See the detail in Appendix of Ref.y]. 


parameters values 

center of wavelength of creation laser Ac 800 nm 

relative line width of creation laser {Suj/ < uj >) 7.5 x 10“^ 

center of wavelength of inducing laser A, 1064 nm 

relative line width of inducing laser < W 4 >) 1.0 x 10 ““^ 

duration time of creation laser pulse per injection Tc 900 fs 

duration time of inducing laser pulse per injection r* 9 ns 

creation laser energy per Tc 9.3 ± 1.2 /rJ 

inducing laser energy per Tj 100 ± 1 /rJ 

focal length / 200 mm 

beam diameter of laser beams d 20 mm 

upper mass range given by 0 < A0 0.15 eV 

u = oji/uj 0.75 

incident-plane-rotation factor Q ^fJ=19/32 

Gl^=l/2 


axially asymmetric factor Tg 


J-ff22=19.4 

-Ff2^12 = 19-2 


combinatorial factor in luminosity Cmb 1/2 

single-photon detection efficiency e^) 1.4 ± 0.1 % 

efficiency of optical path from interaction point to pathjl} Copti 0.5 ± 0.1 % 
efficiency of optical path from interaction point to path{2} eopt 2 0.9 ± 0.2 % 
5Nsi 2.2 

5Ns2 4.4 
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detection efficiency of the PMT mainly caused by the quantum efficiency of the device, eo 
is evaluated using a 532 nm pulse laser in advance of the search. We evaluate the abso¬ 
lute detection efficiency by splitting the 532 nm beam equally and taking the ratio between 
these energies. The one is measured by a calibrated beam energy meter and the other is 
measured by that PMT with neutral density hlters with measured attenuation factors. We 
then corrected the difference of the quantum efficiencies between 532 nm and 641 nm lights 
by taking the relative quantum efficiencies provided by HAMAMATSU into account. 

We then evaluate upper limits on the coupling-mass relation at a 95% conhdence level on 
the basis that the fluctuation of the number of signal yields forms a Gaussian distribution. 
We define 5Ns as the one standard deviation of Ns- It is evaluated from the quadratic sum 
of statistical and systematic errors in Eq. ffT^ and 2.2‘i5Ns is the upper limit of Ns when 


we obtain a 95% conhdence level 


3l[ |. The upper limit of signal yields per shot Wc (for the 


scalar held exchange) and (for the pseudoscalar held exchange) are evaluated as follows: 


3^.c = 


2.2A5Ns2 


yps = 


2.2A5Nsi 


(15) 


^opt2^DWs' ^optl^oWs 

As we briehy mention in Introduction and in detail in Appendix A, even though we hx linear 
polarization planes for creation and inducing laser helds by the polarizers at the moment 
of plane wave propagation, mixing of {1} and {2}-polarization states is unavoidable in 
the focused QPS. By this ehect, the focused system has sensitivity to both scalar and 
pseudoscalar helds simultaneously. 

We obtain the coupling-mass relation from Eq.(IH]). The exclusion limits for scalar and 
pseudoscalar helds at a 95% conhdence level are shown in Fig. [TT]and Fig. [12], respectively. 


VIII. CONCLUSIONS 

A search for scalar and pseudoscalar helds via the four-wave mixing precess at QPS has been 
performed by focusing 10 /iJ/0.9 ps pulse laser and 100 pJ/9 ns pulse lasers. The number of 
{1} and {2}-polarized signal-like photons are Nsi = 0±0(stat.)±2.16(syst.I)±0.30(syst.n)± 
O(syst.ni) and Ns 2 = 1-46 ± 1.27(stat.) ± 2.16(syst.I) ± 0.04(syst.n) ± 3.59(syst.I[l), respec¬ 
tively. We conhrmed that the expected number of four-wave mixing photons in the residual 
gas are negligibly small by measuring the pressure dependence. As a result, no signihcant 
four-wave mixing signal is observed in this experiment. We obtained the upper limits on the 
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FIG. 11: Exclusion limits for scalar fields ((/>) in (/>-photon coupling {g/M) as a function of mass of 
(j) {m^). The excluded region by this experiment is drawn by the red shaded area. The magenta 
shaded area shows the excluded region by our previous search, which is renewed from the black 
dotted line obtained from Ref. 2^ by taking the incident-plane-rotation factor Q and the mass- 


dependent W factor in Appendix B into account. The blue shaded area represents the excluded 
region for scalar fields by light shining through a wall experiment ’’ALPS” 19| (For the mass region 


above 10 ^eV, the sine function part of the sensitivity curve is simplified to unity for drawing 
purposes). The green shaded areas indicate the limits given by non-Newtonian force searches by 


torsion balance experiments ’’Irvine” 


32|. ”Eto-wash 


and Casimir force measurement ’’Lamoreaux” 


],Q, 


’Stanfordl” 


351 ■ ”Stanford2” 


^1 


371]. 


coupling-mass relation for scalar and pseudoscalar fields at a 95% confidence level, respec¬ 
tively. The most sensitive coupling limits g/M = 5.24 x 10“^GeV~^ for scalar search and 
g/M = 5.42 X 10“'^GeV~^ for pseudoscalar search are obtained at m = 0.15 eV. 
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FIG. 12: Exclusion limits for pseudoscalar fields (cr) in cr-photon coupling (g/M) as a function 
of mass of a {ma). The blue shaded area shows the excluded region by the pseudoscalar search, 


’ALPS”. The green and gray solid line show the exc 


’Tokyo Axion Helioscope” 9-11] and ’’CAST 


usion limits from the solar axion experiments 


- IlSl j , respectively. The black shaded area 


represents the result from the dark matter axion search using a microwave cavity ” AD MX” 20|, [2l| . 
The c yan band indicates the expected coupling-mass relation of QCD axion predicted by KSVZ 
model 




39(1 with \E/N — 1.95] in the range 0.07-7, furthermore, in the case of E/N = 0 is shown 


by the black dotted line. 
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Appendix A: Evaluation of the incident-plane-rotation factor Q 



FIG. 13; Definitions of polarization vectors and rotation angles in QPS. 


Figure [13] illustrates the relation between experimentally dehned linear polarization di¬ 
rections {1} and {2} and those theoretically dehned (1) and (2). It also depicts relations 
between pi — p 2 and ps — p^ planes with respect to the x — z plane where the theoretically 
allowed coupling of an exchanged held to the linear polarization states can be evaluated 
in the clearest way. In Ref.j^, we have assumed the incident photons pi and p 2 are both 
plane waves with diherent wave vectors on the same reaction plane which always ensures the 
clearest condition. In the general 3-dimensional incident case such as a focused Gaussian 
beam, however, a pi — p 2 plane can rotate with respect to the x — z plane, which results in 
a deviation from the theoretically clearest condition. We, therefore, introduce a weighted 
averaging factor Q over the clockwise rotation angle <h of the incident reaction plane with 
respect to the x-axis as follows. 

As we have discussed in ref.j^, the Lorentz invariant s-channel scattering amplitude for 
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Lagrangian defined in Eq.Q have the following basic form 




vll'vfJ 


(Pl +V 2 Y + 

where S = abed with a,b,c,d = 1 or 2, respectively, denotes a sequence 
polarization states and m is the mass of scalar or pseudoscalar held. With 
below, the vertex factors for the scalar case (SC) are expressed as 


( 16 ) 

of four-photon 
vectors dehned 


= (PiP2)(eS“^e?^) - (pie^“^)(p2eS'’^), 

= iP3Pi){e^3^ef^) - {P3ef){pief), (17) 


and these for the pseudoscalar case (PS) are expressed as 


^;mps 

^ab 


c PliiP2p^lu ^2a 1 

y[2]PS _ _ pupa 

^cd ~ ^ PSpPApi^Su ^Aa ■ 


(18) 


We must hrst take into account the clockwise rotation angle (p of — Pa plane with respect 
to the given x — z plane independent of the pi — p 2 plane, because these two planes are not 
coplanar in QPS contrary to the situation where the coplanar condition of pi through p^ is 
always satished in CMS. This implies that the simple summation factor 27r on the azimuthal 
degree of freedom of solid angle cannot be applied to QPS, instead, the (p-dependent squared 
transition amplitude must be summed over the possible rotation ip from 0 to 27r. We have 
already introduced this axially asymmetric factor JTs with respect only to the incident 
reaction plane at <h = 0 in j^. This factor essentially depends only on the second vertex 
factors above, while the incident-plane-rotation factor Q is relevant only to the hrst vertex 
factors. We thus dehne the incident-plane-rotation factor as a weighted average with respect 
to J^s at $ = 0 as follows 


/o'|rS(i> = o)N4>’ 


(19) 


because experiments cannot hx the incident reaction plane and intensity of the creation laser 
held must be shared over possible incident reaction planes. 

By requiring (!)={!} and (2)={2} at = (p = 0 where theoretically clearest polarization 
relations can interface with the experimental condition, we describe the polarization vectors 
and momentum vectors for four photons with rotation angles $ and ip as follows: 


eP = (0,1,0), 


( 20 ) 
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= (—COST?, 0, sin??), 62^^ = (—cos??, 0, — sin'i?), 

63^^ = (—cos 03, 0, sin 6*3), 64^^ = (— cos 04, 0, — sin 04), 

Pi = (u;sin'0cos$, —Ci;sin'0sin$,a;cos'0;a;), 

P2 = (—a;sin'0cos$,a;sin'0sin$,a;cos'0;ci;), 

P3 = (a;3sin03cos99,—a;3sin03sin99,ci;3cos03;ci;3), 

P4 = (—a;4sin04cos</9,a;4sin04sin(y9,a;4cos04;a;4). (21) 


We note here that we cannot rotate polarization vectors because the experiment must in¬ 
troduce hxed polarization vectors. This implies that the clear distinction between scalar 
and pseudoscalar couplings cannot be stated due to non-zero rotation angles because non¬ 
identical linear polarization planes between photon 1 and 2 or photon 3 and 4 are implicitly 
introduced. 

Based on these vectors, we summarize relations between momenta and polarization vectors 
with photon labels i = 1, 2, 3,4 as follows 


= —a; slu'd sin $, = a; slu'd sin 4 ), 

= —wa sin 03 sin (p, ^^46^^^^ = a;4sin04sin(p, ( 22 ) 

= 1 and = 0 (23) 

for any pair i, j, and 

=1 for i= j, 

= cos2'd, (63^^64^^^ = cos(03 -|- 04) = cos 0 +, 

= cos('d - 03), = cos('d - 04), 

= cos (^9 -f 04), = cos (^9 -f 03), ( 24 ) 


and 

{P 1 P 2 ) = a;^(cos2'd - 1) = (^ 3 ^ 4 ) = a;3i^4(cos0+ - 1) 
where (pi + P 2 Y = {ps + PaY is required for massless photons. 


(25) 


We are now ready to estimate the factor Q included in the partially integrated cross section 
at Eq.(A24) in Ref. 22|]. We evaluate the case of a6 = 11 for the scalar exchange. From the 


21 



first of Eq. flT7|) . we obtain 


= {piP2){e^i^e^2^) - {pie^2^){p2e^l^) 

= u‘^{cos2'd - 1 + sin^'(9sin^ $) ~ u‘^'d‘^{2 - sin^ $), 


(26) 


where the hrst of Eq. fl2^ . Eq. fEUD . and {pie' 2 ^){p 2 e^i^) = — (wsin'dsind*)^ are substituted. 
The last approximation is based on "d ~'d,, ^ 1- 

This yields the following averaging factor on the incident reaction plane 

/o^’^(2-sin2<h)2d<f> 19 


rsc 
yu = 


'11 


Svr 


32 ■ 


(27) 


We also provide the case of ah = 12 for the pseudoscalar exchange as follows. Based on the 
hrst of Eq. fllSp . the hrst vertex factor with vector dehnitions above is expressed as 


yWPS 

Pl2 


.frMl'PO- p(l)p(2) 

c PltiP2p(ili/ ^2a 


-T)i T)o ci^yp^ 
PlpP2p^ ^2(j 


-PlpP2p 
P2p 


= P2p 


= U 


ePyp-(-cos sin 

{pioe^^^^ + cosi9 + sini9 

(^-coe^yp^ + a; cos cos i? + + a; sin cos sin i9 

(^-uje^y^^P 2 z + uj cos cos i? + + w sin i9 cos sin'd 

^ (—cos'd + cos'd) cos'd+(—slu'd — slu'd) cos $ slu'd = —2a;^ sin^'d cos $. (28) 


This yields the following averaging factor on the incident reaction plane 


r^ps ft cos^ (fdp 

=-2)^- 


I 

2‘ 


(29) 


Appendix B: Refinement of the weight factor W 


In Ref. { 3 , 122 !, we approximated W as a constant 7r/2 for the mass region much smaller than 
that covered by Ad as a conservative estimate. This is because we rather respected simplicity 
of the parametrization than accuracy. However, once we need to compare the sensitivity 
for the higher mass region with the other search methods, the validity of the approximation 
applicable only to the smaller mass region must be reconsidered. In the following, we first 


exactly repeat the relevant part of Ref. 


22 ! and then rehne W as a function of sensitive mass 


regions by quoting necessary equations. 
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We first express the squared scattering amplitude for the case when a low-mass field is 
exchanged in the s-channel via a resonance state with the symbol to describe polarization 
combinations of initial and hnal states S. 


^ you 

X + 0, 

where y = with the resonance condition m = 2ur sin for a given mass m and a 


is expressed as 


/gm\‘^ 


Stt V M 


mT 

2 sin^ -dr 


with the resonance decay rate of the low-mass held 


r = (167r)-^ 


The resonance condition is satished when the center-of-mass system (CMS) energy between 
incident two photons Eqms = 2 c<;sin-d coincides with the given mass m. At a focused 
geometry of an incident laser beam, however, Eqms cannot be uniquely specihed due to 
the momentum uncertainty of incident waves. Although the incident laser energy has the 
intrinsic uncertainty, the momentnm uncertainty or the angnlar uncertainty between a pair 
of incident photons dominates that of the incident energy. Therefore, we consider the case 
where only angles of incidence -d between randomly chosen pairs of photons are nncertain 
within 0 < -d < A-d for a given focusing parameter by hxing the incident energy. The 
treatment for the intrinsic energy uncertainty is explained in Appendix B later. We £x the 
laser energy u at the optical wavelength 


< = ~ (33) 

while the resonance condition depends on the incident angle nncertainty. This gives the 
expression for y as a fnnction of -d 


/ON 2 2/ ON 

X(^) = = 


(l - Ujlpt, 


where 


dd = 


X X- 




We thus introduce the averaging process for the squared amplitude IAT 5 P over the possible 
uncertainty on incident angles 

\M^= r^%(d)|Af 5 (^)Pdd (36) 
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where Ms specified with a set of physical parameters m and gM~^ is expressed as a function 
of'd, and p('d) is the probability distribution function as a function of the uncertainty on 
within an incident pulse. 

We review the expression for the electric field of the Gaussian laser propagating along the 


^-direction in spatial coordinates {x, y, z) 
^(x,y,z) = Eq 


as follows: 


Wq 


exp < —i[kz — H{z)] — 


+ 


ik 


(37) 


where Eq electric field amplitude, k = r = + y‘^, wq is the minimum waist, which 

cannot be smaller than A due to the diffraction limit, and other definitions are as follows: 


w{zf=wo^l + ^], 




(38) 


R = z\l + 


zr^ 


(39) 


H{z) = tan ( — ) , 
zr 


(40) 


Zr = 


TTtCo 


A 


(41) 


With 9 being an incident angle of a single photon in the Gaussian beam, the angular 
distribution g{9) can be approximated as 

1 1 


9{0) 


exp 


02 I 

‘2A02 j ’ 


(42) 


V^Ae 

where the incident angle uncertainty in the Gaussian beam A6 is introduced within the 
physical range |0| < 12 as 


A0 


A. 


d 


(43) 


ITMIo 2/’ 

with the wavelength of the creation laser Ac, the beam diameter d, the focal length /, and 
the beam waist wq = as illustrated in FigUl For a pair of photons 1, 2 each of which 
follows g{9), the incident angle between them is defined as 


^ = ^\ei-e2\. 


(44) 
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With the variance A’d'^ = 2{jA9‘^), the pair angular distribution p('d) is then approximated 


as 


pW 


y/lTAd 


exp < — 



y/irAd 


for 0 < ■(? < 7r/2 


(45) 


where the coefficient 2 of the amplitude is caused by limiting d to the range 0 < "d < 7r/2, 
and -C 1 is taken into account because A9 in Eq. (l4^ also corresponds to the upper 

limit by the focusing lens based on geometric optics. This distribution is consistent with the 
flat top distribution applied to Ref.jsl, [2^ except the coefficient. 

We now re-express the average of the squared scattering amplitude as a function of y = 
in units of the width of the Breit-Wigner(BW) distribution a by substituting Eq. dSU]) and 
into Eq. fl36|) with Eq.( 


f 1 



[ a 9 ) 


where we introduce the following constant 


aW, 


(46) 


>V = 






e + 1 




(47) 


with 

W(0 = (1 - (48) 

^opt 

In Eq. fHT)) the weight function W{^) is the positive and monotonic function within the 
integral range and the second term is the Breit-Wigner(BW) function with the width of 
unity. Note that is now explicitly proportional to a but not a?. This gives the 

enhancement factor a compared to the case jWlsP oc where no resonance state is contained 
in the integral range controlled by A9 experimentally. The integrated value of the pure BW 
function from ^ = —1 to ^ = -|-1 gives 7r/2, while that from ^ = —oo to ^ = -|-oo gives vr. 
The difference is only a factor of two. The weight function 1E(.^) of the kernel is almost 
unity for small a^, that is, when a is small enough with a small mass and a weak coupling. 
Therefore, we will consider only the region of ± 1 as a conservative estimate. By taking 
only this integral range, we can be released from trivial numerical modihcations originating 
from ^ = —oo and the behavior of W(^) at ^ = -^{1 — {'dr/ (7r/2))^} which are not essential 
due to the strong suppression by the Breit-Wigner weight. 

We now rehne W in order to apply it more accurately even to the case for Pr/A9 ~ 1 
where, exactly speaking, the second approximation in Eq. fl45|) is not valid. In this case, by 
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using the first of Eq. (H5|) with substitution of the relation between x = and'd expressed 
in Eq. fl3i|) . Eq. fHHjl is modihed as follows 



(49) 


where the last approximation is based on <C 1 with respect to the integral range ^ ± 1 


in Eq. fl47p for the conservative estimate. This is justihed in the mass-coupling range we are 
interested in via the hrst relation in Eq. flST]) . for instance, ~ 10“^® for m ~ 0.1 eV 

and g/M GeV“^. By substituting Eq. fH^ into Eq. fH7|) . the conservative evaluation 

on W over ± 1 is expressed as 



(50) 


This factor is dependent ofequivalently dependent of mass, especially for larger close 
to A6 while it is almost 7i/2 for smaller 
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